Theoretical energies and transition probabilities of argon 
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Argon energies of J=l odd and J=2 even states, oscillator strengths of the two lowest J=l odd 
states, and line strengths for transitions between J=l odd and J=2 even excited states are calculated 
with the CI+MBPT method. The results are in good agreement with experiment. 
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I. INTRODUCTION 

Given that a combined configuration-interaction (CI) 
and many-body-perturbation-theory (MBPT) approach 
could describe well experimental data for energies and 
transition rates in neon, it was a natural continuation 
to investigate a more difficult, stronger core-correlated 
atom such as argon. In argon theories and experiments 
appear to be less accurate than in neon, with larger mu- 
tual disagreement; moreover, heavier closed-shell atoms 
are even less studied and understood so that the argon 
atom can serve as an intermediate step in understanding 
many closed-shell atoms and ions. 

Although semiempirical calculations Q, 0, 0, 0> 
sometimes lead to a good agreement with experiment for 
some transition arrays and provide valuable information 
for the analysis of existing experimental data, they do 
not explain much the nature of correlations and their 
predictive power is limited to cases where (a) energies 
are well-known and (b) energy-level information is suffi- 
cient to infer other atomic properties such as transition 
probabilities. In general two theories could give exactly 
the same energy levels without predicting the same tran- 
sition probabilities. For example, adding or not adding 
random-phase approximation (RPA) corrections to tran- 
sition amplitudes will give two different answers for tran- 
sition rates although wave functions obtained from the 
same effective Hamiltonian are identical. When RPA 
corrections were included in neon |fj, the improvement 
of agreement with experiment was only marginal and it 
was not an issue, but as we will show in this paper this 
type correction is crucial in argon. Hence, theories, which 
do not take into account at least the RPA correction, can 
be misleading in argon and heavier noble gas atoms. Still 
owing to cancellations of RPA contributions for the same 
transition array, ratios obtained with a matrix method 
could be still accurate. Indeed, a fair agreement with 
experiment for low-lying argon transitions is obtained 
by Lilly Q even without considering RPA corrections. 
However, the accuracy of ratios decreases if relativistic 
and correlation effects, that introduce asymmetry for the 
P3/2 an d P1/2 hole RPA contributions to radial transition 



integrals, are significant. 

Another issue is the role of configuration interaction 
since many theories are limited by considering only a few 
configurations. The convergence of CI has to be investi- 
gated. Furthermore, other many-body effects, such as an 
interaction of a valence electron with a core or a hole with 
a core or their combination, are quite significant in argon 
and heavier noble-gas atoms. Thus, the development of 
ab initio theory for heavy noble-gas atoms is needed to 
single out important correlation and relativistic effects. 
Some success is already achieved in this direction. For ex- 
ample, elaborate ab initio calculations by Avgoustoglou 
and Beck performed for transitions from the lowest 



odd excited states 



3 P3/2 4s l/2 



and 



3 Pi/2 4s i/2 



to the 
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ground state give oscillator strengths in good agreement 
with average experimental values. However, this theory 
was not extended to transitions between excited states. 
Therefore, one objective of this paper is to prove that 
transitions between excited states of argon can also be 
calculated with ab initio methods. We also would like 
to show that for argon our simpler CI+MBPT theory 
has comparable accuracy to complicated single-double 
couple-cluster (SDCC) method of Ref. Q and thus can 
be of practical value for other theoretical groups. 

In this paper, first we will briefly describe our method 
of calculations (more details are given in H, @); then, 
we will calculate energies and oscillator strengths of J=l 
odd states and energies of 3=2 even states and compare 
them with experiment. This comparison gives an esti- 
mate on accuracy of our wave functions. Next, we will 
present line strengths and rates for transitions between 
excited states and compare our result with experiment 
and semiempirical calculations. 



II. CI+MBPT METHOD 

We have developed a convergent variant of CI+MBPT 
for particle-hole states of closed-shell atoms [t| and suc- 
cessfully applied it in neon to calculate energies § 
and oscillator strengths for transitions between excited 
states 0- This perturbation theory is based on the first 
iteration of couple-cluster single-double equations [Tol 
ITU . All second-order MBPT terms [l^] are included, 
with some denominators being modified to take into ac- 
count the strong interaction between a hole and a core 



2 



electron or core and core electrons non-perurbatively |8|, 
9]. The advantage of this approach compared to the it- 
erated single-double method 11] is simplicity and speed 
of calculations. The accuracy of hole energies and fine- 
structure splittings is significantly improved already af- 
ter adding second-order MBPT corrections as was illus- 
trated previously in neon. Apart from Coulomb correla- 
tion corrections, the Breit magnetic interaction is also in- 
cluded, but small frequency-dependent Breit, quantum- 
electrodynamic, reduced-mass, and mass-polarization 
corrections are omitted. 

To calculate particle-hole energies, we construct a 
model CI space, compute effective Hamiltonian in this 
space, and solve an eigenvalue problem. Along with en- 
ergies we obtain wave functions, which are used to cal- 
culate transition amplitudes or other related quantities. 
We investigate energies, oscillator strengths for transi- 
tions to the ground state, and line strengths or transition 
rates for transitions between excited states. Energies of 
neon particle- hole J=l odd states and oscillator strengths 
were in very good agreement with experiment: pure ab- 
initio neon energies differed from experimental energies 
by 0.0069 a.u., but after subtraction of the systematic 
shift (which does not make much difference in transition 
calculations), the agreement was at the level of 0.0001 
a.u. for almost all states. In this paper we use the same 
approach to calculate energies and transition amplitudes 
of argon. 

Reduced matrix elements for transitions between 
particle- hole excited states are given in Ref. |(J: 



(F \\Zj\\ I) = y/(2J F + 1) (2 J/ + 1) 



|(-i) j + j '+^+jv 



J I J Jp 

Jv' Ja Jv 



Sa'a (v'\\Zj\\v) + 



Jv Ja' Ja 



We found that second-order corrections to transition am- 
plitudes or lowest order RPA corrections 



7RPA 



E 

bm 



J brn Qimjb 



E 

bm 



Qibjr 



(1) 



TABLE I: Energies of argon odd J=l states. For unique spec- 
ification of levels, accurate NIST energies are provided. The 
sizes of CI model spaces are denoted CI-22, etc. A7 is the 
deviation from experiment for the 7th level. All energies are 
expressed in atomic units 



State NIST CI-22 CI-32 CI-52 



A 7 



P 3 /2 4S 



Pi/ 

P 3 /2 3d 

Pgy 2 5S 

P 3 ~/2 3d 
Pl/2 5s 

pr/ 2 3d 



0.4272 
0.4347 
0.5095 
0.5178 
0.5201 
0.5239 
0.5256 



0.3202 
0.4219 
0.4295 
0.5082 
0.5234 
0.5268 
0.5331 



0.4219 
0.4291 
0.4348 
0.5085 
0.5257 
0.5268 
0.5330 



0.4219 
0.4297 
0.5083 
0.5250 
0.5268 
0.5317 
0.5332 



0.0053 
0.0050 
0.0012 
-0.0072 
-0.0067 
-0.0078 
-0.0076 



0.0129 
0.0126 
0.0088 
0.0004 
0.0009 
-0.0002 
0.0000 



denominators is important and brings our theory to a 
good agreement with experiment which will be illustrated 
later in this paper. 



III. CALCULATIONS FOR ARGON 

A. Energies of argon odd J=l states 

Table [Q shows our calculated energies for J=l odd 
states. Fine structure for the 4s, 5s, or 4f groups of states 
is well reproduced since the deviation from experiment 
for each fine structure group of levels is similar. The last 
column shows shifted deviations from experiment, which 
demonstrate that there is a systematic shift for energy 
levels in our calculations. This shift is due to the inac- 
curacy of our calculations of hole energies. The shifted 
deviations illustrate a good precision for higher excited 
states which have smaller correlation effects. We also 
show the convergence with the number of particle-hole 
configurations. The change in energy is significant when 
we increase the size of CI space from 22 to 32 and to 
52; however, from 52 to 62 the change is small, and we 
consider that CI-52 case is optimal. 



are significant in argon. This type correction requires 
careful analysis since the possible interaction between an 
"observing" hole state a and core states b over which 
summation is carried out should be taken into account in 
all orders, at least a bulk part of this interaction. Similar 
to our approach for energies, we include the interaction 
in all orders by adding the extra term 



extra 



Xo(abab) 



y/(2j a + l)(2j b + l) 



(2) 



in the denominators of RPA expressions. This procedure 
reduces RPA corrections roughly twice. If we note that 
the value of RPA correction is comparable to the values 
of lowest order matrix elements, the modification of the 



B. Oscillator strengths of resonant transitions 

We calculate oscillator strengths for well-studied ex- 
perimentally and theoretically resonant argon transi- 
tions. The upper and lower panels of Fig. ^ show data 
for the [3^3/2 4s] 1 and [3p]y 2 4s]i states, respectively. The 
agreement with SD CC calculations and most experi- 
ments is very close. For the weaker transition, deviation 
from the SD CC theory is somewhat larger, which is due 
to cancellation of contributions from different configu- 
rations. Such cancellation was more complete in neon, 
where we observed larger disagreement with the theory 
of Ref. 0. 
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TABLE II: References for experimental data shown in Fig. 



Experiments 

CI+MBPT 

SDCC 



T 

T---J- 



T 



6 8 10 12 

Observations 
Ar[3p 1 1/2 4s]i 



14 



16 



18 



Experiments 
- CI+MBPT 
SDCC 



6 8 10 12 

Observations 



14 



16 



18 



FIG. 1: Comparison of CI+MBPT with experiment and SD 
CC calculations Q for oscillator strengths of the [p^^sji and 

[pT^sJi states of argon 



C. Energies of argon even J=2 states 

Table 11111 shows our calculated energies for J=2 even 
states. Results for these and J=l odd states are similar: 
fine structure for 4p, 5p, or 4f states is well reproduced; a 
small systematic shift exists due hole energy offset which 
if subtracted brings theory in a better agreement for less 
correlated higher excited states. The convergence with 
the number of particle-hole configurations is achieved for 
the size of CI space 48. 



D. Argon transitions between excited states 



Table llVl shows the comparison of our calculations with 
experiment for line strengths between J=l odd and J=2 
even state transitions. Very large disagreement is ob- 
served when no RPA correction is added. However, if 
conventional RPA correction (Eq^l is included, the dif- 
ference is also large, though of an opposite sign. And only 



# Reference 



Year 



/ Is* 



/ lS2 



1 Lewis [13] 

2 Lawrence ^4] 

3 Geiger [15] 

4 de Jongh and Eck [16] 

5 Irwin et al. [17] 

6 McConkey and 
Donaldson 
Copley and Camm 
Vallee et al. [20] 
Westerveld et al. [21] 

10 Hahn and 
Schwentner [22J 

11 Chornay et al. [23] 

12 Li et al. [24] 

13 Tsurubuchi et al. [25] 

14 Chan et al. [26] 

15 Ligtenberg et al. [27] 

14 Wu et al. [28] 

15 Gibson and Risley [29] 



1967 
1968 
1970 
1971 

1973 0.083±0.027 



0.063±0.003 
0.059±0.003 
0.047±0.009 



7 



9 



1973 
1974 
1977 
1979 

1980 
1984 
1988 
1990 
1992 
1994 
1995 
1995 



0.096±0.002 
0.076±0.008 
0.051±0.007 
0.063±0.005 



0.077±0. 

0.065±0. 

0.0581^ 

0.057±o! 
0.0662±0. 
0.0616±0. 
0.0676±0. 
0.0580±0. 



004 
005 

005 
008 

003 

0033 

0021 

0040 

0017 



0.278±0.020 
0.228±0.021 
0.186±0.037 
0.22±0.02 
0.35±0.13 



0.283±0.024 
0.210±0.030 
0.240±0.020 

0.361±0.060 

n 999 +0.002 
U.z^z_ 003 

0.213±0.011 
0.265±0.013 

0.2297±0.0093 
0.259±0.015 

0.2214±0.( 



TABLE III: Energies of argon even J=2 states. For unique 
specification of levels, accurate NIST energies are provided. 
All energies are expressed in atomic units. A7 is the deviation 
from experiment for the 7th level 



State 

P3/2 4 P 
P3/2 4 P 

Pi~/2 4 P 

P3/2 5 P 
P3/2 5 P 
Pl"/2 5 P 

?3/2' 



NIST CI-6 CI-24 CI- 



A — A 



7 



P,",n4f 



P3/2 4f 
P 3 /2 6 P 



0.48123 
0.48405 
0.48885 
0.53309 
0.53393 
0.53978 
0.54762 
0.54781 
0.55219 



0.4838 
0.4870 
0.4919 
0.5397 
0.5405 
0.5466 



0.4837 
0.4871 
0.4920 
0.5396 
0.5406 
0.5466 
0.5597 
0.5601 
0.5665 



0.4844 
0.4878 
0.4928 
0.5403 
0.5413 
0.5474 
0.5562 
0.5564 
0.5604 



0.0032 
0.0037 
0.0039 
0.0072 
0.0074 
0.0076 
0.0086 
0.0086 
0.0082 



-0.0050 
-0.0045 
-0.0043 
-0.0010 
-0.0008 
-0.0006 
0.0004 
0.0004 
0.0000 



if we include the appropriate, modified-denominator RPA 
correction (mRPA), which constitutes roughly a half of 
the conventional RPA correction, do our values agree 
with experiment at the level of the experimental accu- 
racy. The experimental accuracy in Table llVl is denoted 
by letters in accordance with the NIST convention: for 
example, the accuracy of class "C" is about 10%. Still in 
the case of the p^^s-p^^p transition the disagreement 
is quite large and independent of RPA corrections which 
might indicate that the problem is in the experiment. If 
we exclude this case and the inaccurate measurement of 
class "D" , then mRPA calculations can be considered 
accurate at the level of 10-20%. 

We also compare our final results (mRPA calculations) 
with some available semiempirical calculations by Lilly 



4 



TABLE IV: Line strengths in atomic units for transitions 
between J=l odd and J=2 even states. Accuracy of ex- 
periments is shown as letters in accordance with the NIST 
convention [30||. Three cases with respect to RPA correc- 
tions are considered: "cRPA" (conventional), "nRPA" (no), 
and "mRPA" (modified-denominator) RPA corrections are in- 
cluded. The deviation from experiment is measured by A, 
which is defined as A = (S cxp . - Stheory )/S e x P . 



5.0 -, 



4.5 



4.0- 



4s- 
4s- 



4s- 



'3/2 L 
'3/2 E 



'3/2 ' 



1/2 



„4s- 



1/2 



'3/2 L 



Pl/2- 
Pl/2 4S -P372 f 



LI 



LL 



8427 


31.8 


C 


24.7 


31.6 


27.8 


13% 


8008 


6.21 


C 


2.60 


3.28 


2.94 


53% 


7386 


8.42 


C 


7.31 


9.41 


8.23 


2% 


4301 


0.074 


C 


0.016 


0.15 


0.062 


-16% 


4267 


0.060 


C 


0.025 


0.079 


0.048 


20% 


4046 


0.054 


C 


0.025 


0.099 


0.054 


0% 


9787 


3.40 


c 


3.41 


4.26 


3.79 


-12% 


9227 


9.75 


c 


10.9 


13.8 


12.2 


-25% 


8411 


32.7 


c 


22.7 


28.9 


25.6 


22% 


4630 0.0094 C- 


0.0013 


0.024 


0.0097 


-3% 


4591 


0.0015 


D 


0.0003 


0.0041 


0.0006 


-60% 


4335 


0.114 


c 


0.030 


0.214 


0.0985 


14% 



3.5 



Transition A,A Exp. Ac. cRPA nRPA mRPA A~ < 
pTpS^pp 8497 31 8 n 247 3T~fi 97 8 13% 3.0 

„-l A„ „-l A„ 
^3/2 

P3/2 

P3/2 

P3/2 

Pi/2 



2.5 



2.0 



1 Experiments 
— Current theory 
NIST 



-T--I- 



l 



T 



4 6 8 10 

Observations in chronological order 



i — 
12 



14 



FIG. 2: Comparison of CI-MBPT and experimental transition 
probabilities for the 430-nm argon transition 



TABLE VI: References for experimental data shown in Fig. [21 



TABLE V: Comparison of transition rates for transitions be- 
tween J=l odd and J=2 even states. NIST values represent 
experimental results, our transition rates are obtained from 
mRPA line strenths in the previous table using experimental 
wavelengths. Wavelengths specify transitions shown in the 
previous table. In addition to experiment, our results are 
also compared with various semiempirical calculations. Rates 
are expressed in s _1 . Brackets denote powers of 10 



A, A NIST mRPA Ref.fl 



8427 2.15[7 

8008 4.90[6 : 

7386 8.47[6 

4301 3.77[5 : 

4267 3.13[5 : 

4046 3.30[5 

9787 1.47[6 : 

9227 5.03[6 

8411 2.23 [7 

4630 3.84[4 

4591 6.28[3 

4335 5.67[5 



1.88 [7" 
2.32[6 

8. 28 [6 
3. 16 [5 
2. 50 [5 
3.30[5 
1.64[6 

6. 29 [6 
1.74 [7" 
3.96[4 
2.51[3 
4.90[5 



2.08[7 
4.93[6 
9.18[6 : 
1.33 [6 
3.90[5 
5.30[5 
1.25[6 : 
5.24[6 : 
2.25 [7 
2.20[5 
1.50[5 : 
1.66[6 : 



Ref.[2 



2.11[7 

4.8 [6 ; 
9.11[6 ; 

1.3[6 : 
4.6[5 
5.2[5 
1.07[6 : 

5.9 [6 ; 
2.21 [7 

2.0[5 : 
1.6[5 
1.7[6 



Ref.[3] Ref.[4] Ref.[5] 



3. 77 [5] 
2.80[5] 
3.40[5] 



3.9[4] 
3. 60 [3] 
5. 90 [6] 



3.77[5] 3.77[5] 
3.51[5] 3.05[5] 
3.10[5] 3.49[5] 



3.7[4] 3.9[4] 
7[3] 8[3] 
5.53[6] 5.62[6] 



H, Garstang and Blerkom 0, Wende 0], Bues et al. 
U, Jones and Wiese 1. Income cases semiempirical 
results are close to experiment, but in others they can 
deviate substantially, by orders of magnitude, and they 
also disagree with each other. For example, for weak 
transitions, transition rates of Lilly [| are about 10 times 
larger than experimental values; for these transitions re- 
sults of Garstang and Blerkom [!| are much better. Our 
theory agrees with experiment more consistently in cases 



Obs. 


Reference 


Year A (10 b s _i ) 


a 


1 


Drawin [31] 


1956 


3.58 


20% 


2 


Gericke [3J 


1961 


3.25 


15% 


3 


Richter [33] 


1965 


3.1 


10% 


4 


Popenoe and Shumaker, Jr. \34j 


1965 


4.11 


11% 


5 


Chapelle et al. [3^ 


1968 


3.14 


12% 


6 


Wende [3J 


1968 


3.1 


18% 


7 


Wujec [3J 


1971 


3.66 


10% 


8 


van Houwelingen and Kruithof [3JJ 


1971 


3.08 


8.5% 


9 


Nubbemeyer [3^ 


1976 


3.40 


5% 


10 


Preston [39] 


1977 


3.72 


4% 


11 


Baessler and Kock [40] 


1980 


3.2 


7% 


12 


Hirabayashi et al. [4L] 


1987 


3.94 


15% 



where experiments have good precision. 

For one particular transition at 430 nm, many mea- 
surements have been conducted, and this transition can 
be considered to be known most accurately among other 
transitions between excited states of argon. In Fig. [21 we 
compare our theoretical value with many available ex- 
perimental values. While six experiments support NIST 
quoted value, a greater number of experiments support 
our theoretical value. Therefore, we tend to concluded 
that NIST values can be somewhat wrong. We also 
find that semiempirical calculations are not able to pre- 
dict these particular experiments; namely, most exten- 
sive semiempirical calculations |l| give the transition rate 
1.33 x 10 6 s _1 that disagrees with all experiments. Re- 
turning to our Table ITVl we can argue that some large 
deviations of our theory from experiment can be caused 
by inaccuracy of quoted experimental errors. 

For better tests of our theory, more accurate experi- 



5 



merits are needed. For this reason, we have measured 
the ratio of line strengths using an accurate laser absorp- 
tion technique Our experimental ratio, 3.29 ± 0.13, 
for the two transitions at 9787 A and 9227 A completely 
agrees with our theoretical value, 3.22. It would be ex- 
tremely interesting if this experimental technique was 
further developed and applied to other transitions in ar- 
gon or in other noble gas atoms. For example, experi- 
ments in neon using visible lasers could test the accuracy 
of this experimental method, since our theory is expected 
to be more accurate in neon than in argon. 



particle-hole states. Our theoretical energies and transi- 
tion rates are in good agreement with most experiments. 
We found that random-phase approximation corrections 
are essential. We also analyzed dependence of accuracy 
on the size of configuration space and determined the op- 
timal size of CI matrix. Our CI+MBPT program can be 
used for calculations of many other transitions of closed- 
shell atoms and ions. 



IV. CONCLUSIONS 



In this paper, we have applied CI+MBPT theory 
to calculate energies and transition properties of argon 
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